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Abstract. We investigate geometric aspects of the the Backlund trans- 
form of principal contact element nets. A Backlund transform exists if 

— and only if it the principal contact element net is of constant negative 

Gaussian curvature (a pseudosphere). We describe an elementary con- 
f--) struction of the Backlund transform and prove its correctness. Finally, we 

show that Bianchi's Permutation Theorem remains valid in our discrete 

, * setting. 
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1. Introduction 

r— i In the 1880s A. V. Backlund and L. Bianchi explored the transformation 

of one surface of constant negative Gaussian curvature (a pseudosphere) 
Q into a surface of the same constant Gaussian curvature. These surface 

transformations were named after Backlund and, in their formulation in the 
language of partial differential equations, play an important role in soliton 
theory and in integrable systems. We refer the reader to the monograph 
i— i [ 10 1 for a comprehensive modern treatment with many applications. 

In 1952, W. Wunderlich gave a geometric description of the Backlund 
^ transform of a discrete structure that is nowadays called a K-net — a discrete 

ON asymptotic net of constant negative Gaussian curvature |IT4| . Later, analytic 

formulations were added and the discrete transformations were embedded 
ff~j into the theory of Discrete Differential Geometry, see the monograph (4|. 

q In this article we perform a comprehensive geometric study of the Back- 

lund transform of principal contact element nets — nets of contact elements 
such that any two neighbouring contact elements have a common tangent 
sphere. Our main results are original but there exist related contributions 
. ^ in the above mentioned publications. We also have to mention the article 

^ |[TT| by W. Schief. It contains a description of the Backlund transform of 

discrete O-surfaces. While the vertex sets of discrete O-surfaces and princi- 
pal contact element nets are identical (both are circular nets), the normals 
differ. For O-surfaces, they are defined by a simple algebraic condition, for 
principal contact element nets they satisfy a geometric criterion. Accord- 
ingly, our approach is of geometric nature as opposed to Schief 's analytic 
treatment. 

Our results are natural extensions of recent works in discrete kinematics 
[[121 [13| . While their formulation is rather natural at this point, some of 
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our proofs tend to be rather involved and require extensive calculations. 
We partly attribute this to the observation that a (smooth or discrete) 
asymptotic parametrization is more appropriate for the description of 
Backlund transforms. Principal contact element nets discretize principal 
parametrizations and lead to spatial structures of complicated nature. 

The occasional use of a computer algebra system (CAS) is indispens- 
able in this work. We use Maple 13 for this purpose. In order to make 
the computer calculations understandable, we try to give rather explicit 
descriptions. 

2. Preliminaries and statement of main results 

We proceed by giving definitions for the main concepts, by stating our 
central results and by clarifying their relations. Proofs are deferred to later 
sections. 

Principal contact element nets provide a rich discrete surface representa- 
tion, consisting of points and normal vectors. They have been introduced 
in [3] in an attempt to develop a common master theory for circular nets 
(see [4, Section 3.5]) and conical nets (see [4, Section 3.4] and ll9lP. 

Definition 1. A contact element is a pair (p, n) consisting of a point a G E 3 
(Euclidean three-space) and a unit vector n 6 M. 3 . Its normal line N is the 
oriented straight line through a and in direction of n, its oriented tangent 
plane n is the plane through p with oriented normal vector n. 

In this text we will always denote normal vector and normal line by 
the same letter but in different fonts: Lowercase, boldface for vectors, 
uppercase, italic for lines — just as in Definition [T] 

Definition 2. A contact element net is a map (p,n) from Z d , d ^ 2, to the 
space of contact elements. A principal contact element net is a contact element 
net such that any two neighbouring contact elements (pi,ni) and (pj,nj) 
have a common oriented tangent sphere. 

The defining condition of principal contact element nets is rather re- 
strictive. It implies that neighbouring normal lines N| and Nj intersect 
in a common point h|j which is at the same oriented distance from both 
vertices pt and pj . This shows that neighbouring contact elements have a 
bisector plane (3 ij . Moreover, the vertices of an elementary quadrilateral lie 
on a circle, the tangent planes are tangent to a cone of revolution and the 
normal lines form a skew quadrilateral on a hyperboloid of revolution. 

An elementary quadrilateral (po,no), (pi,ni), (p2,n2), (p3,ns) of a prin- 
cipal contact element net can be constructed by choosing the vertices po, 
Pi/ Vi> and P3 on a circle, prescribin no and then finding the other normal 
vectors by reflection in the bisector planes |3i,i+i of pt and pi+i- 

In this article we investigate principal contact element nets that admit a 
Backlund transform: 

Definition 3. Two principal contact element nets (p,n) and (q, m) are called 
Backlund mates if 
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(1) The distance of corresponding points pi and qi is constant. 

(2) The angle of corresponding normals and mt is constant. (It is 
necessary to measure this angle consistently in counter-clockwise 
direction when viewed along the ray from pi to qt.) 

(3) Two corresponding tangent planes 7tt and Xi intersect in the con- 
necting line pi V qi of their vertices. 

In this case, (q,m) is also called a Backlund transform of (p,n) (see Figure [TJ. 

The last condition in Definition [3] could be replaced by the requirement 
that the connecting line pi V qi is perpendicular to both normal vectors, ni 
and rrii. 

In the classical theory, Backlund transforms of smooth surfaces O, W 
can be defined by the same three conditions as in Definition [3] (with the 
understanding that both surfaces are parametrized over the same domain 
and corresponding points, normals, and tangent planes belong to the same 
parameter values). It is our aim in this article to prove several results on 
Backlund transforms of principal contact element nets that are also valid 
in the smooth theory and in other discrete settings. The most important 
property states that only pseudospheres admit Backlund transforms. 

Theorem 4. If a principal contact element net (p,n) admits a Backlund transform, 
it is of constant negative Gaussian curvature 

m k = -^ 

where d is the distance between corresponding points and cp is the angle between 
corresponding normals. 

We still have to explain the notion of Gaussian curvature to which this 
theorem refers. In the smooth setting, the Gaussian curvature in a point 
can be defined as the local area distortion under the Gauss map. This is 
imitated in 
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Figure 2. Transversals of four relative revolute axes 

Definition 5. The Gaussian curvature K of an elementary quadrilateral 
(po,no), (pi,ni), (p2,n2), (p3,n3) of a principal contact element net is 
defined as 

(2) K = § 

where So is the oriented area of the spherical quadrilateral no, ni, n2, n3 
and S is the oriented area of the circular quadrilateral po, pi, P2/ P3- 

The Gaussian curvature as defined here is a well-accepted concept in 
discrete differential geometry (see [4]) and it is a special case of a recently 
published and more general theory 0. Note that in contrast to the smooth 
setting, the Gaussian curvature is assigned to a face and not to a vertex. This 
allows us to speak of the Gaussian curvature of an elementary quadrilateral. 

It will be convenient to speak of the twist of two oriented lines and also 
of the twist of Backlund mates. 

Definition 6. The twist of two oriented lines N and M is defined as the 
ratio 

(3) twist(N,M) = 

d 

where (p is the angle between the respective direction vectors n and m 
(measured according to the convention of Definition [3jl and d is the distance 
between N and D (Figure [2j. The twist of two Backlund mates is defined as 
the twist of any two corresponding normal lines. 

We will argue that Theorem [4] is a consequence of the Theorem |HJ below, 
which deserves interest in its own. In order to formulate this result, we 
recall a definition from [13|: 

Definition 7. A cyclic sequence Aq, Ai, A2, A3 of direct Euclidean displace- 
ments Ai : E 3 — > E 3 is called a rotation quadrilateral, if any two neighboring 
positions correspond in a rotation. 

A generic rotation quadrilateral has four relative revolute axes Kqi, K12, 
K23, K30 (in the moving space) that admit two transversal lines N and M. 
Denote the feet of their common perpendicular by p G N and q € M and 
one of their unit direction vectors by n and m (thus imposing an orientation 
on N and M). It was shown in Ifl3| that precisely the homologous images 
of the two contact elements (p,n) and (q, m) can serve as corresponding 
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elementary quadrilaterals of Backlund mates. This leads us to an important 
result that relates the Backlund transform of pseudospherical principal 
contact element nets to discrete rotating motions. 

Theorem 8. Consider a rotation quadrilateral Aq, A\, A2, A3 with relative 
revolute axes Kqi, K12, K23, K30 in the moving space. Denote by N and M their 
transversals, by n and m one of their respective unit direction vectors and by 
p £ N and q G M the feet of their common normals (Figure^. Then the Gaussian 
curvature of the homologous images of (p, n) and (q, m) equals — twist 2 (N, M), 
that is, it can be computed by ^ where d = dist(p, q) and (p = <(n,m). 

As a consequence of Theorem [8] a pseudospherical principal contact 
element net of Gaussian curvature K admits Backlund transforms of twist 




In retrospect and especially when studying its relation to Theorem |4j the 
content of Theorem [B] is maybe a rather obvious conjecture. Nonetheless, 
it is surprising when viewed independently. There are six degrees of 
freedom when fitting a rotation quadrilateral to two given lines N and M: 
Three choices for every relative rotation makes a total of twelve degrees 
of freedom; the requirement that their composition yields the identity 
consumes six of them. Still, this is insufficient to change the area ratio 

Our proof of Theorem [B] uses the dual quaternion calculus of rigid mo- 
tions. It allows a computational proof that requires only linear operations. 
Nonetheless, the involved calculations are so complicated that assistance of 
a CAS seems indispensable. The reason for this is the rather involved spa- 
tial configuration and the high number of degrees of freedom. Theorems |4] 
and [FJ will both be proved in Section 3] 

A certain converse of Theorem |4j is the construction of the Backlund 
transform from a given pseudospherical principal contact element net. It is 
on the agenda in Section [4] 

Theorem 9. Given a pseudospherical principal contact element net (p,n) of 
negative Gaussian curvature K and a contact element (q^mO such that 

• qi lies in the tangent plane 7tt o/pt, 

• mi is perpendicular to the line pt V qt, and 

• the quantities K, d = dist(pt, qt) and cp = <(nt,mt) satisfy ([I) 

there exists precisely one Backlund transform (q,m) o/(p,n) such that (pt,ru) 
and (qt,mt) correspond. 

As we will see, there exists a simple necessary construction for the 
neighbours of the contact element (qt, mt). This already implies uniqueness 
of the Backlund transform. In order to prove existence, we have to show 
that the conditions on qt and mt imply that this construction does not 
lead to contradictions. The proof consists once more of a direct CAS-aided 
computation. 

One interpretation of Theorem [9] is that every pseudospherical principal 
contact element net can be generated in infinitely many ways as trajectory 
of a discrete rotating motion that has a second, non-parallel trajectory 
surface of the same type. Combining this observation with results of 
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fl2| , in particular with |l2l Theorem 7], we get an interesting corollary 
(the terminology is explained in the standard reference [4] on discrete 
differential geometry): 

Corollary 10. Principal contact element nets of a prescribed negative Gaussian 
curvature are described by a multidimensionally consistent ID-system. 

This corollary certainly admits simple direct proofs. Still, we are not 
aware of a reference that actually states it. 

We should also mention that our results are in perfect analogy to the 
smooth setting. If <D and W are Backlund mates, the correspondence 
between their points can be realized by principal parametrizations x(u,v) 
and y(u,v). For every parameter pair (u,v), we obtain a rigid figure 
consisting of corresponding points x(u, v), y (u, v) and the respective surface 
normals. This induces a two-parametric motion A(u,v) with the following 
properties: 

• A(u,v) is a gliding motion on both surfaces O and W [8, Sec- 
tion 7.1.5]. This means that there exist two planes in the moving 
space whose images under A(u,v) envelop O and W, respectively. 

• The infinitesimal motions in u- and v-parameter directions are 
infinitesimal rotations. (This is a general property of gliding motions 
along principally parametrized surfaces.) 

In our discrete setting, both properties are preserved. The infinitesimal 
rotations of the continuous case are replaced by rotations between finitely 
separated positions. 

Our final result, which will be presented in more detail in Section |5j is 
Bianchi's Permutation Theorem for Backlund transforms of pseudospheri- 
cal principal contact element nets: 

Theorem 11. Suppose (b,l) and (c,m) are two Backlund transforms of the same 
twist to a principal contact element net (a,k). Then there exists a unique principal 
contact element net (d, n) which is at the same time a Backlund transform of (b, 1) 
and (c,m). 

Note that throughout this article we implicitely make a number of 
assumptions on the genericity of the involved geometric entities. More 
specifically, we require that the Gaussian curvature Q is well-defined for 
all elementary quadrilaterals of principal contact element nets (with the 
notable exception of Example [15| T he meaning of the word "generic" in 
the paragraph foiling Definition U\ is not as easily explained. We require 
that indeed two transversals to the four relative revolute axes exist in 
algebraic sense. This can be formulated as a non-vanishing condition of 
an algebraic expression in the input paramters. Moreover, neighbouring 
relative revolute axes are assumed to be skew. 

3. The Gaussian curvature of Backlund mates 

This section is dedicated to the proof of Theorem |4] Assume that (p, n) 
and (q, m) are Backlund mates, denote the distance of corresponding points 
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by d and the angle of corresponding normals by cp. We want to show that 
the Gaussian curvature of all elementary quadrilaterals of (p,n) is constant 
and can be computed by Q. 

Our first observation concerns two neighbouring pairs of corresponding 
contact elements. 

Lemma 12. Under the conditions of Theorem |£j any two neighbouring pairs 
(po,rto), (pi,rti) and (qo,rrto), (qi,rrti) of corresponding contact elements corre- 
spond in a rotation Rqi about the line Lqi = hoi V kpi where hoi = Nq n Ni and 
koi = M n Mi. 

Proof. By Definition |3j the figures formed by corresponding contact el- 
ements are congruent and of equal orientation. Thus, there exists a 
direct Euclidean displacement Roi such that Roi(po,no) = (pi,ni) and 
Roi(qo/mo) = (qi,mi). By Definition [2] the points hoi and koi are fix-points 
of Roi- Thus, Roi is indeed a rotation about the line Loi. □ 

Consider now corresponding elementary quadrilaterals (po,no), 
(p3,n3), and (qo,mo), . . . , (q3,rri3) of the two Backhand mates. They define 
four relative rotations Roi, Ri2/ R23, R30 such that 

(4) R 30 o R 23 o Ri 2 o R i = id. 

This observation relates the Backlund transform of pseudospherical princi- 
pal contact element nets to the geometry of rotation quadrilaterals |13| and 
discrete gliding motions fll"2) . We will use results and techniques from both 
articles for our proof. 

More specifically, we view each figure (pi,ni), (qi,rru) as the position 
of a rigid body in space. This gives us four direct Euclidean displacements 

(5) A| : E 3 — )- E 3 

from the moving space E 3 to an identical copy of E 3 , the fixed space. They 
form a rotation quadrilateral in the sense of Definition [7] Theorem [8] states 
that the Gaussian curvature of an elementary quadrilateral depends only 
on the relative position of two corresponding normal lines. Since this is the 
same for all elementary quadrilaterals, Theorem [4] follows. 

In our proof of Theorem [8] we use the dual quaternion calculus of spatial 
kinematics. It is explained for example in 0. A short primer is also 
given in [12]. We use the notation of the latter article. A direct Euclidean 
displacement is modelled as a unit dual quaternion A = A + eA where 
A and A are its respective primal and dual part and e satisifes e 2 = 0. 
Assembling the components of primal and dual parts into a homogeneous 
vector [A, A] = [ao, . . . , ay], we can view A as a point on the 

(6) ScP 7 :A-A = 0. 
This quadric is called the Study quadric. 

Proof of Theorem |8] We start by setting up a system of equations that de- 
scribe a rotation quadrilateral whose relative rotation axes in the moving 
space intersect two fixed lines N and M. At first, we specify the two lines 
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N and M in the moving space. Without loss of generality, the normal feet 
p and q on N and M, respectively, can be chosen as 

(7) p = (0,0,e), q = (0,0,-e) 

where 2e = d. The respective directions vectors of N and M are 

(8) n=(l,t,0), m = (l,-t,0) 

where t is related to the angle cp between N and M via cp = 2 arctan t. Note 
that the vectors n and m are not normalized so that we ultimately will not 
compute the Gaussian image of an elementary quadrilateral on the unit 
sphere but on a sphere of squared radius 1 + 1 2 . This is admissible since the 
Gaussian curvature will only be multiplied by the constant factor 1 + 1 2 . 

Without loss of generality, the first position of the rotation quadrilat- 
eral can be taken as the identity, Ao = [1, 0, ... , 0] (we use homogeneous 
coordinates in P 7 ). The two neighbouring positions are 

(9) Ai = [ai ,...,ai 7 ], and A 3 = [a 30 , . . . , a 37 ]. 
Their components ciij are subjects to seven constraints: 

• The relative displacements Roi and R03 are rotations. In the dual 
quaternion calculus this constraint is modelled as 

(10) n 5 {A i *A 1 )=0, i€{l,3}, 

where 715 denotes the projection onto the fifth coordinate. 

• The relative revolute axes K01 and K03 intersect M and N: 

(11) a((Roi) e *B(Roi*A ) e *(R^ 1 ) £/ T)=0; i e {1,3}, T e {M, N}. 

In this equation, D. is the bilinear form associated to the Study 
quadric (Q(X, Y) = X • Y), the subscript e denotes e -conjugation of 
dual quaternions, and t)X is the vector part of the dual quaternion X. 

• The sought positions lie on the Study quadric §: 

(12) Q(A, A) = Ai • At = 0, i G {1, 3}. 

Equations ( |T0"| l and (IT) are linear in the components of Ai and A3, Equa- 
tion ( |T2"} is quadratic. Normalizing by aio = tai3, the general solution for 
Ai reads 

^ Ai = (tai3,taiiai3,ta 2 3 ,0,-t 2 anai3, Qi2ai3,t 2 an - a 2 2 ); 

ail/ an, ai3 G R. 

Here, an, a.12, and 013 serve as free parameters. The solution for A3 is 
obtained by replacing the parameter aij with a^j (j = 1,2,3). 

Similarly, the missing position A2 is uniquely determined by two linear 
equations Ei, E2 of type ( fTO) , four linear equation E 3 , . . . , Eg of type ( flT) and 
one quadratic equation E7 of type dl2j. Thus, it is the second intersection 
point of a straight line through Aq with the Study quadric S. The solution 
is unique and can be computed in rational arithmetic. 

While it is possible to compute A2 and verify ([l) directly, the involved 
expressions are excessively long. Our implementation of this approach 
delivers the desired result but only after a few hours of computation. 
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Therefore, we favor a method which is based on ideal theory and yields 
a confirmatory answer within a few minutesj^] The basic idea is to show 
that a certain polynomial (Equation ( [181 , below) is contained in the ideal 
spanned by the polynomial equations Ei, . . . , E7. 

For i G {0, 1,2,3} we compute the images pt of the normal foot p and the 
images nt of the direction vector n under the displacement At (with the 
entries a.2i of A2 left unspecified). This can be accomplished by using the 
following formulas for the action of a dual quaternion A on a homogeneous 
coordinate vector [xo, X\, %2, X3] = [xo,x]: 

(14) x' Q + ex' = A e * (x + ex) * A 

(the bar denotes dual quaternion conjugation). Note that the use of ho- 
mogeneous coordinates avoids the introduction of denominators — this is 
highly desirable in CAS calculations of a large complexity — but has to be 
taken into account later in ( [15} when computing oriented areas. 

Instead of computing the oriented area of the quadrilaterals po, pi, p2/ 
P3 and no, ni, n2, n3 we compute the areas S and So of their respective 
projection onto the first and second coordinate. This is possible since an 
orthographic projection does not affect the ratio of areas in parallel planes. 
We have 

3 

S = Y (PUPi+1,2 -Pi / 2Pi+l / l) x i+2><i+3/ 
i=0 
3 

So = y (nj,l n i+l,2 — n i / 2rtj+l / l) x i+2^ r i+3 



(15) 



(indices modulo four). The factors x\ = YS)=o a ij compensate the use of 
homogeneous coordinates in the computation of the points pi and the ideal 
points (vectors) n^. Now we have to show that 



(16) 



Sn t 2 



S (l + t 2 )e 2 ' 
or, equivalently 

(17) (l + t 2 )e 2 S -t 2 S =0. 



Equation 17 is polynomial. It holds true if its left-hand side is contained 
in the ideal J spanned by the seven equations Ei, . . . , E7 that determine 
the fourth position of A2. The indeterminates are CI20, ■ ■ ■ , 0.27 while an, 
ii2/ CL13, a.31, 032, 033, e, t serve as parameters. We perform the necessary 
algebraic manipulations by means of the CAS Maple 13 and assume that the 
defining equations of A2 are stored in El, . . . , E7. Moreover, the simplifying 
normalization (120 = 1 is applied. 
with(Groebner) : 

# J defined by linear and quadratic equations 
J := [El, E2, E3, E4, E5, E6, E7] : 



We thank Dominic Walter for his help with this approach. 
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# Use tdeg term order: 

TO := tdeg(a21, a22, a23, a24, a25, a26, a27) : 

# Interreduce J and reduce R with respect to this reduced ideal. 

# The result is indeed 0, showing that R is contained in J. 
J := InterReduce( J, TO): 

Reduce ( [R] , J, TO) [1] ; 

The final output is 0, showing that R is indeed contained in the ideal J. □ 

This proof also finishes the proof of Theorem [4] A few explaining 
remarks concerning the computer calculation seem appropriate: 

• The tdeg-term ordering (Line 6) is primarily by total degree and 
then reversed lexicographic. 

• Inter-reducing the generating elements of the ideal J with respect 
to the term order TO (Line 10) produces a list of polynomials that 
generate the same ideal J but are reduced in the sense that no 
monomial is reducible by the leading monomial of another list 
element. 

• The command Reduce (Line 11) computes the remainder of R di- 
vided by the polynomials in J. Its vanishing implies that R is indeed 
contained in the ideal spanned by J. 

The running time of this implementation is approximately six minutes. 
Minor improvements are possibly by adapting the order of the variables in 
Line 6. We regret not being able to offer a proof of Theorem [4] that does not 
rely on a computer algebra system. The main reason for the computational 
difficulties (in spite of the fact that the problem is linear) is the high number 
of six free parameters. A more insightful proof or a proof with manually 
tractable calculations would be desirable. 

4. Construction of the Backlund transform 

In this section we provide an algorithm for actually constructing (or 
computing) the Backlund transforms of a pseudospherical principal contact 
element net. The proof of correctness of this algorithm is at the same time 
the proof of Theorem [9] 

It is a simple observation that the Backlund transform (q,m) in The- 
orem [9] is necessarily unique. Existence is a different matter. Given the 
principal contact element net (p, n) and the contact element (q^mi), any 
neighbour (q^m-j) to (q^mi) is uniquely determined. It is found by the 
following steps whose necessity has already been discussed: 

• Determine the points h.^ = N| n Nj and k|j = n |3ij where (3ij 
is the bisector plane of pi and pj. 

• Determine the unique rotation R|j about the axis Ltj = Ktj V ktj that 
transforms (pt,nt) into (pj,nj). 

• The contact element (qj,rrij) is the image of (q-^mt) under the 
rotation . 
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We refer to this construction briefly as the "neighbour construction". Note 
that a minor variation takes as input the oriented lines N-t, Nj and Mi (as 
opposed to three contact elements) and returns the oriented line Mj . 

We provide some more information on the actual computation of the 
neighbour construction — not only as a service for the reader but also 
because we need them later in the proof of Lemma 13 and Theorem [9] 

Throughout this calculation we use homogeneous coordinates to describe 
points and planes, Pliicker coordinates to describe lines and homogeneous 
four by four matrices to describe Euclidean displacements. (Here, the 
use of dual quaternion calculus does not provide a significant advantage.) 
The bisector plane (3 of two non-ideal points [xo, X\, X2, X3] = [xo,x] and 
[yo/"y 1/^2/^3] = tyo/y] nas homogeneous coordinates 

(18) [uoX-x-Xoy-y,2x -yoxxy]. 

The intersection point of the plane [tio, Ui, U2, U3] = [uo, u] and the line with 
Pliicker coordinates [gi,g2/93/hi, h.2, h.3] = [g,h] equals 

(19) [x , x] = [u • g, -uog + u x h] . 

Because of hij = Nt n |3ij = Nj n (3^ this formula can be used to determine 
both, k|j and hy . 

The rotation R^j can be conveniently written as the composition of two 
reflections in the bisector plane (3ij and the plane y^ = pj V hy Vkij. 
This is true since the composition of two reflections in planes |3 and y is 
a rotation about the intersection line L = (3 n y. Moreover, it is obvious 
that (pi, n-t) is transformed into (pj,nj). The homogeneous transformation 
matrix [Rij] of Rij is given by 

(20) [R tj ] = [C tj ] • [By] 

where [Bij] and [Cij] are the reflection matrices in (3ij and yij, respectively. 
They can be computed by means of the general formula for the reflection 
matrix in a plane [uq, vl\, U2, U3] which reads 

"uf + u^ + u^ 

— 2uoUi — u^ + u| + u| — 2uiU2 — 2U1U3 

— 2uqU2 — 2U}U2 Uj — + u| — 2U2U3 

— 2U0U3 — 2U1U3 — 2U2U3 + — Ug_ 

Equations (|l8"|l-(|2"T|| are sufficient for a CAS implementation of the neigh- 
bour construction. 

We return to the proof of Theorem [9] where existence is still open. The 
problem with the neighbour construction is that it will produce contra- 
dictions for a general choice of (qo, mo). Consider only an elementary 
quadrilateral (po,no), (p3,n3) and the contact element (qo,mo). By 
the neighbour construction we obtain (in that order) the contact elements 
(qi,mi), (q2,rri2), and (q3,rri3). Applying the neighbour construction once 
more, we get a contact element (q^rrig) which should equal (qo,mo). 

Since the Gaussian curvature K and the distance d between correspond- 
ing points is already defined, we necessarily have to choose mo so that 
the angle cp = <(no,mo) satisfies ([lj. Thus, there are only four choices 



(21) 
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Figure 3. Notation in Lemma [13] 

for the unit vector mo and only two choices for the normal line Mo- We 
will show that all of them lead to valid solutions. This is also enough to 
ensure that the neighbour construction can be consistently applied to the 
whole principal contact element net. The proof of this still needs some 
preparatory work. 

To begin with, we mention those configuration that will be identified as 
false positives by our test (qo,mo) = (q^m^). 

• If Mo intersects the axis of the circle through po, pi, P2, P3, the com- 
position of the four rotations yields the identity but the homologous 
contact elements (qo,mo), . . . , ^3,1113) do not form the elementary 
quadrilateral of a principal contact element net. 

• If Mo is parallel to No the action of the rotation R 1,1+1 on (qi,mi) 
equals that of the reflection in Therefore, we get (q^m^) = 
(qo,mo) also in this case. The Gaussian curvature of the two ele- 
mentary quadrilaterals is, however, different. 

It goes without saying that both configurations violate the pre-requisites of 
Theorem [HJ 

The following lemma states that the neighbour construction acts projec- 
tively on the set of lines obtained by revolving Mo about No- It allows a 
simplifying assumption during the later calculations. 

Lemma 13. Consider three lines Nq, Ni and Mq such that Mq and Nq are skew 
while Nq and Ni are intersecting. Denote the set of lines obtained by rotating 
Mq about Nq by Mo- Then the neighbour construction with respect to Nq and Ni 
(equipped with some orientation) induces a projective mapping between Mo and 
its image Mi under the neighbour construction. (Figure^. 

Proof. Our computation is based on a Cartesian coordinate frame {o;x,y,z} 
which is projectively extended such that vanishing of the first coordinate 
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characterizes ideal points. We choose No as the z-axis and Mo such that 
the shortest distance between No and Mo is attained at the pedal points 
[1,0,0,0] G N and [1,1,0,0] € M . The ideal point on M than equals 
[0, 0, 1, v] with some v£l. Because of the rotational symmetry of the set 
Mo, we may assign the coordinates [s,0,t,u] to the bisector plane (3oi of Nq 
and N\. Identifying geometric entities with their homogeneous coordinate 
vectors, we can write 

N = [0,0,1,0,0,0], Ni = [0,2tu,u 2 -t 2 ,2st,0,0], 

(22) hoi = No n |3oi = No n Ni = [u,0,0,-s], 

Mo(A) = [-2A,A 2 -l,u(A 2 + l),2uA,u(l-A 2 ),l + A 2 ]. 

As expected, Mo (A) is a quadratic parametrization in Pliicker coordinates 
of the line-set Mo- The theorem's statement amounts to saying that Mi (A) 
is also quadratic. The intersection point of Mo with (3oi equals 

uv-t+ (t + uv)A 2 
t - uv + 2sA + (t + uv)A 2 
s + 2uvA — sA 2 
— sv — 2tvA — svA 2 



(23) koi(A) 



Now we compute Mi (A) by reflecting Mo (A), at first in the plane (3oi 
and then in yoi = Ni Vkoi(A). We refrain from giving the details of 



the calculation; the relevant formulas are ( 18 ) — <[2Tj>. The outcome is a 
polynomial vector function Mi (A) which a priory is of degree 20 in A. 
However, all of its entries have the common factor 

(24) (t 2 + u 2 ) 9 (l + A 2 ) 8 G(A) 
with the quadratic polynomial 

(25) G (A) = s 2 + 1 2 - 2tuv + uV + 4stA + (s 2 + 1 2 + 2tuv + u 2 v 2 )A 2 . 



This common factor ( 24 1 can be split off so that only a quadratic component 



remains. □ 



Remark 14. The vanishing of the polynomial ( 24 1 has a geometric meaning. 

• If t 2 + u 2 vanishes, the plane (3oi is at infinity and the reflection in 
|3oi becomes undefined. 

• If 1 + A 2 vanishes, the direction of the revolute axis (3oi n Yoi is 
isotropic and the rotation becomes singular. The same is true for 
the zeros of G(A). 

All of these instances are irrelevant in our setting. 

Lemma [13] implies that the map 
(26) rj : M -»• M , M i-> MJ, 

as a composition of neighbour maps, is projective as well. The line-set Mq 
is of course the set obtained by subjecting Mq to all rotations about Nq. 
We conclude that r| generically has two fixed lines — those lines of Mo that 
intersect the z-axis. We are interested in configurations where a non-trivial 
fixed line Mq G Mq exists. In this case n is the identity on M. This can also 
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be phrased as follows: If (qo,mo) is such that the neighbour construction is free 
of contradictions, the same is true for every contact element that is obtained from 
(qo,mo) by a rotation about No. This observation gives us one extra degree 
of freedom for a simplifying assumption in the 

Proof of Theorem [9] The basic idea is to show that in the pencil of lines Mo 
incident with qo and perpendicular to po V qo precisely two lines result in 
a closed chain of neighbouring contact elements obtained by successive 
application of the neighbour construction. These solutions necessarily 
satisfy the condition of the theorem. 

We begin by assigning homogeneous coordinates to the four points po, 
Pi/ P2/ P3- Without loss of generality we may set 



(27) 



Pi = [l + til-tf^O], i = 0,... ,3. 



The values t| are in E U {oo} and ti = oo corresponds to the point [1, — 1, 0, 0]. 
The following calculations do not take into account this exceptional value 
but could easily be adapted to handle this situation. Moreover, we make 
the admissible simplification to = 0. 

The bisector planes of the points pt and pj are 



(28) 



where = t-t + tj and T|j 
the matrix 



(3ij = [0,5^,^,0], t,j =0,...,3; i^j 

titj — 1. The reflection in (3ij is described by 



(29) 



[Bt 



S?. + T- 2 . 

^ij 'ij 











-I? 



s 2 . 





^StjTtj 







s 2 . 











S 2 . + T 2 . 
^ij + 'ij. 



Starting with the normal vector no = [0,u,v, 1], we compute 
(30) ni = [B i] • n , n 3 = [B 30 ] • n , n 2 = [B 2 i] • n x . 

The first point qo of the Backhand transform should be such that the 



connecting line po V qo is perpendicular to no. Due to Lemma 13 the 
possible choices for mo can be parametrized as 



(31) 



mo = [0, Au, Av, ev + A] . 



Note that the invalid choice mo || no corresponds to A — > oo and will 
automatically drop out during our calculations. The second invalid choice 
is obtained for A = in which case the line Mo intersects the circle axis. 
The point qo is then 



(32) 



q = [1 + 1§, 1 - 1§ + ev,2t - eu,0] 



where e is a parameter that determines the distance between po and qo- 
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Using the additional admissible simplification to = 0, the intersection 
point koj = Mo n (3oi becomes 

(uti — v)A 
-(v + e(u 2 +v 2 ))A 
-(v + e(u 2 + v 2 ))tiA 
-(ti +tie(u + v))(ev + A) 



(33) k 0j 



j €{1,3}. 



By means of ( |20| and pTj ), we compute the rotations Roi and R03 to obtain 
the homogeneous coordinates of the points qi, q 3 and the vectors (ideal 
points) mi, 1113. Common factors of the shape 

(34) eV+v^l+tJ) 4 , j€{l,3}, 

can be split off. The resulting expressions are just a little too long to be 
displayed here. 

In the same way, we compute the intersection points k.12 and k23, where 
we immediately split off the common factor 

(35) (l + t 2 ) 2 ((u 2 +v 2 + l)(utj -v) 2 A 2 +v 2 (eu + t j (ev + l)) 2 ) / )e{l,3}, 

the rotation R12, and the vector rri2. This latter vector can be divided by 

e 2 (u 2 +v 2 )(l + t 2 ) 6 (l + t 2 ) 4 Q(A) where 

(36) Q(A) = (u 2 +v 2 + l) 2 (tiu-v) 4 A 4 + 

2v 2 (eu + t!(ev + l)) 2 (u 2 + v 2 + l)(tiu-v) 2 A 2 +v 4 (eu + ti(ev + l)) 4 . 

The closure condition of the neighbor construction is now the linear 
dependence of rri2 and the ideal point (direction vector) [0, k 2 ] of k.12 V IC23. 
It turns out that k 2 can be divided by 

(37) (t 1 -t 3 )(v + e(u 2 + v 2 ))A. 

We compute the cross-product rri2 x ]<2 and extract the greatest common 
divisor 

(38) 2evt 2 Gi(A)G 2 (A) 

of its entries. Here, Gi(A) and G2(A) are quadratic polynomials in A. 

One of them, say Gi, equals A -1 times the homogenizing coordinate of 
ki2- It is a spurious solution since it vanishes in configurations where 1<2 
and rri2 are parallel but M2 is different from ki2 V k23. The second factor 
equals 

(39) G 2 =v 2 (T3(e 2 u 2 + l-e 2 v 2 )-2T 2 e 2 uv) + +2ev(T3(u 2 -v 2 )-2T 2 uv)A 

+ (T 3 (u 2 - v 2 ) - 2uvT 2 ) (v 2 + u 2 + 1)A 2 

where 

(40) T 2 = t 1 t2-tit3 + t 2 t3 + l and T 3 = tit 2 t 3 + t t - 1 2 + 1 3 . 

Its vanishing characterizes the positions of Mo that lead to valid solutions. 
Since |39| is quadratic in A, there are precisely two solutions. This finishes 
the proof of Theorem [9] □ 
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Example 15. The maybe simplest example of a pair of Backhand mates in 
the continuous settings consists of a straight line Z, viewed as a degenerate 
pseudosphere, and the surface of revolution m obtained by rotating a 
tractrix with asymptote Z about Z. We apply our construction to a discrete 
version of this configuration which is actually too irregular to be covered 
by our theory. Nonetheless, we are able to illustrate our basic ideas and to 
recover the geometric essence of the continuous case. 

We define the first principal contact element net (p, n) as follows: 

p:Z 2 ^R 3 , (i,j) ^ (0,0,)) 

(41) 

n: 1} -> R 3 , (i,j) ' y (0,cos(^),sin(^f),0) 

The integer k £ N is a shape parameter that affects the discretization of the 
revolute surface. 

The principal contact element net (p, n) is irregular in the sense that the 
Gaussian curvature of the elementary quadrilaterals is undefined. Nonethe- 
less, we can construct the Backhand transform (q,m) to (p, n), defined by 
the contact element (qo,0/ m o,o) with qo,o = (0, d, 0), mo,o = (0, cos a, sin a) 
and d G M. It is depicted in Figure [4] where we use a = j. 

According to our theory, the contact element (qo,j+i,moj+i) is obtained 
from (p ,j,n ,j), (p ,j + i,nQ,j+i) and (q 0/ j,m 0/ j), see Figure |t} left (where we 
write pj, qj etc. instead of poj, qo,j)- It never leaves the plane x = and 
produces a discrete curve with vertices qoj and normal vectors mo,j. The 
choice a = j ensures that the length of the tangent segment between qoj 
and the z-axis is constant. Hence, we may address the curve as discrete 
tractrix. 

The contact element (qt+i^mi+ij) is constructed from the contact ele- 
ments (py/iiij), (pi+i^ru+y), and (q^my). The intersection point of 
My with the bisector plane of py and Vi+ij is the ideal point of the z-axis 
so that (qi ; j,mt j) and (qi+i^mt+ij) correspond in a rotation about the 
z-axis through the angle The resulting discrete revolute surface can 
be seen in Figure |4j right. Its Gaussian curvature is indeed constant and 
negative. We refrain from deriving analytic expressions for describing this 
discrete surface. These have already been established in HJ. 



5. Bianchi's Permutation Theorem 



Now we prove Bianchi's Permutation Theorem 11 our final major result 



We split our proof into a series of intermediate steps. 

Lemma 16. Assume that (b,l), and (c,m) are both Backlund transforms of the 
same twist o/(a,k). If the contact elements (bo,lo] and (co,mo) correspond to 
(ao,lo) there exists a half-turn that interchanges (bo,lo) and (co,mo) (Figure^. 

Proof. In a suitable Euclidean coordinate frame, we have 

lo = (cos (p, sin (p, 0), mo = (cos cp, — sin cp, 0), 

(42) 

b = (0,0,z) + (31 , c = (0,0,-z)+ym 
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Figure 4. Discrete tractrix and pseudospherical principal 
contact element net of revolution 




Figure 5. Opposite normals correspond in a half-turn 



with (p g {0,j) and (3, y ^ 0. We try to recover the position of (ao/ko). It 
will turn out that this is only possible, if |3 = y, that is, (t>o,lo) and (co,mo) 
correspond in a half-turn. 

Clearly, the possible location of do is the intersection line of the two 
tangent planes to (bo,lo) and (co,mo). It can be parametrized as 

(43) ao(t) = — ((|3+y)sincp,(|3-y)coscp,t). 

2 cos cp sin (p 

The squared distances form ao(t) to bo and Co equal 



(44) 



4 



(Oo-Cq) • (Ofl-Co)/ 



respectively. (We omit the argument t for sake of readability.) Given do, the 
unit normal vector ko is obtained as ko 



| kg 1 1 where 



(45) 



(a -b ) x (d -c ). 
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The squared sines of the angles between the contact element normals are 

(46) sin 2 i|j b = ||1 x k || 2 , sin 2 ^ c = ||m x k || 2 , 

A necessary condition for ao(t) is now the equality of the ratios 

(d7 s sin 2 ih _ sin 2 \L> 3 

' d 2 d 2 

a i Q 3 

or vanishing of the numerator P of 

(48) d 2 sin 2 ih-d 2 sin 2 i|j3. 

Generically P is a polynomial of degree four in t. It can be factored as 

(49) P=4(|3 2 -y 2 )tan 2 |p 1 (t)P 2 (t) 

with quadratic polynomials Pi(t) and P2U) whose explicit form can be 
readily computed using only rational arithmetic. We are rather interested 
in their discriminants. Using the substitution cp = 2arctanu we obtain 

discrim(Pi,t) = -64(1 -u 2 ) 2 u 2 ((|3 -y)(u 4 + 1) + 2(3y + |3)u 2 ) 2 , 

(50) 

discrim(P 3/ t) = -64(1 -u 2 ) 2 u 2 ((|3 -y)(u 4 + 1) - 2(3(3 +y)u 2 ) . 

Thus, no real solutions exist for ao(t) unless (3 2 = y 2 . This must indeed 
be the case, since (a,k) is a real solution. Moreover, since (3 and y are not 
negative, we have (3 = y and the half-turn about the first coordinate axes 
indeed interchanges (b, m) and (c, m). □ 

Lemma 17. Consider three contact elements (ao,ko), (bo,lo), and (co,mo) such 
that twist(Ao, Bo) = twist(Ao, Co). Then there exists a unique contact element 
(do,n ) such that 

• dist(b ,d ) =dist(c ,d ) =dist(a ,b ) = dist(a , c ), 

• no is perpendicular to the vectors that connect do to both, bo and Co, and 

• twist(B , D ) = twist(C , D ) = twist(A , B ) = twist(A , C ). 

Proof. For the time being, we ignore the twist conditions. By the same cal- 



culation as in the proof of Lemma 16 (but with (3 = y) we get distance and 



angle conditions that determine the positions of the contact element (do, no). 
In fact, the possible locus of points do is a straight line, parametrized by a 



vector function do (t) as in (43 1. The unit normal vector no(t) can be found 
as in ( j45) . For all contact elements (do(t),rio(t)) the twist with respect to 
(bo,lo) and (co,mo) is the same. 

Thus, only one twist condition, say twist(Bo, Do) = twist(Ao, Bo), is 
relevant. Its square turns out to be quadratic in t. Thus, there is exactly 
one solution to twist(Bo, Do) = ±twist(Ao, Bo) apart from (ao,ko). It is 
indeed a solution with equal sign, since it can be obtained by rotating 
(ao,ko) through 180° about the half-turn axis that interchanges (bo,lo) and 
(c ,m ). □ 



Lemma 17 allows an unambiguous contact element wise construction of 
the contact element net (d,n). It remains to be shown that (d,n) is indeed 
a principal contact element net. Our proof makes use of a remarkable 
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relation between the Backhand transform and the Bennett linkage that 
was already noticed by Wunderlich in |14|. The original references to the 
Bennett linkage are HJ|2|. A more accessible description of its geometry is 
H Section 10.5]. 

The Bennett linkage is a spatial four-bar mechanism with a one-parametric 
mobility. It consists of four skew revolute axes Ko, Lq, No, Mq such that the 
normal feet on one axes to the two neighbouring axes coincide. Denote 
these points by ao, bo, do, and Cq, respectively. In order to be mobile, the 
mechanism has to fulfill the constraints 

• dist(a ,b ) = dist(c ,d ), dist(a ,c ) = dist(b ,d ), 

• twist(A ,B ) =twist(Co,D ),twist(A ,Co) =twist(B ,D ) 

(see for example [6. Section 10.5]). As suggested by our notation, corre- 



sponding normal lines in Theorem 11 form the axes of a special Bennett 



linkage where all four twists are not only equal in pairs but equal as a 
whole. Thus, we may refer the reader to Figure [5] for an illustration of a 
Bennett linkage. 

A Bennett linkage allows infinitely many incongruent realizations in 
space that exhibit the same relative positions, characterized by normal 
distance and angle, between any two adjacent axes. If one link is kept fix, 
the configuration space of the opposite link has the topology the projective 
line (it can be described as a conic on the Study quadric, see [7|). Hence, 
any two realizations of Bennett's linkage can be continuously transformed 
into each other without changing the relative position of neighboring axes. 
This observation is a key ingredient in the proof of Theorem 



11 



As the final preparatory step for the proof of Theorem 1 1 we mention 
an obvious alternative to the characterization of principal contact element 
nets in Definition [2] Observe that two contact elements (ao,ko) and (bo,lo) 
in general position correspond in a unique rotation. The rotation axis is 
found as line of intersection of two bisector planes of corresponding points 
(for example of ao, bo and ao + ko, bo + lo)- Two rotations exist if and 
only if the bisector planes of two corresponding points coincide. In this 
case the bisector planes of all corresponding points coincide and (ao,ko) 
can be rotated in infinitely many ways to (bo, lo). The rotation axes lie in 
the common bisector plane of corresponding points and are incident with 
Ko n Lo. This leads to 

Proposition 18. A contact element net (a, k) is a principal contact element net if 
and only if any two neighbouring contact elements (at, ki) and (aj, kj ) correspond 
in two (and hence in infinitely many) rotations. The rotation axes are incident with 
the intersection point K| n Kj of the two normals and lie in the unique bisector 
plane of (a^kO and (a^kj). 

Proof of Theorem |TTj We consider an elementary quadrilateral (ao, ko), . . . , 
(a3,k3) of the principal contact element net (a, k) and the corresponding 
quadrilaterals (b ,lo), (b 3 ,l 3 ), and (c ,m ), (c 3 ,m 3 ) under the 



two Backlund transforms. By Lemma 17 there exist uniquely determined 



contact elements (do, no), (d 3 ,n 3 ) that satisfy all distance and angle 
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Ra / 



"0 



lo' 



CO 



mo 




Figure 6. Rotations about Ra (top) and about Li (bottom) 



in the proof of Theorem 11 



constraints imposed by the Backlund transform. We have to show that these 
contact elements form the elementary quadrilateral of a principal contact 
element net. In view of Proposition [18] it is sufficient to show that one 
neighbouring pair, say (do, no) and (di,ni), corresponds in two different 
rotations. 

As already argued earlier, there exists a rotation R that maps (ao,ko) to 
(ai,ki) and, at the same time, (bo,lo) to (bi,li). Its axis Ar is spanned by 
K n Ki and L n L x (Figure [6} top). 

The rotation R transforms (co, mp) to a contact element (c^,m^) and 
(do, no) to a contact element (d^,n^). Since Ki, I_i, Nq, Mq and Ki, Li, Ni, 
Mi are just two realizations of the axes of the same Bennett linkage (with 
respective normal feet ai, bi, dp, Cq and ai, bi, di, Ci), we can rotate (d*, n^) 
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about Ki into (di,ni). (Via Bennett's linkage, this rotation can be coupled 
with a rotation of Cg to C\ about the axis Ki — but this is not relevant at this 
point of our proof.) The composition of R with this last rotation (Figure |6j 
bottom), call it Q, is again a rotation S because the axes Ra and I4 of R 
and Q intersect. Thus, we have found one rotation that transforms (do, no) 
to (di,ni). If we start with the rotation that maps (ao,ko) to (ai,ki) and 
(co,mo) to (ci,mi) we obtain a second rotation T of (do, no) to (di,ni). 
Generically, the two rotations S and T are different, because Lo and Mo do 
not intersect. We conclude that the contact elements (do, no) and (di,ni) 



satisfy the principal contact element criterion of Proposition 18 □ 



6. Conclusion 

In this article we gave a fairly complete geometric treatment of the 
Backlund transform of principal contact element nets. We proved the 
most relevant results (Theorems [5J |9j and [TTj l which are typical of this 
curious surface relation. Moreover, we consider the disclosed relations to 
discrete kinematics and in particular to discrete rotating motions to be of 
interest. In this context we would emphasize Theorem |8j which already 
brings us to open issues. We were unable to provide certain proofs (those 
of Theorems [8] and |9]l and in a way that make the involved calculations 
manually tractable. There is still room left for simplifications. Moreover, an 
analytic complement to our geometric reasoning, maybe in the style of |lTT| . 
would be desirable. 
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